Phys. Rev. B, Rap. Comm. (to appear) 
Decoupling in the ID frustrated quantum XY model and 
Josephson junction ladders: Ising critical behavior 



Enzo Granato 
Laboratorio Associado de Sensores e Materials, 
Institute Nacional de Pesquisas Espaciais, 
12225 Sao Jose dos Campos, S.P. Brazil 

Abstract 

A generalization of the one-dimensional frustrated quantum XY model is 
considered in which the inter and intra-chain coupling constants of the two 
infinite XY (planar rotor) chains have different strengths. The model can de- 
scribe the superconductor to insulator transition due to charging effects in a 
ladder of Josephson junctions in a magnetic field with half a flux quantum per 
plaquette. From a fluctuation-effective action, this transition is expected to 
be in the universality class of the two-dimensional classical XY-Ising model. 
The critical behavior is studied using a Monte Carlo transfer matrix applied 
to the path-integral representation of the model and a finite-size-scaling anal- 
ysis of data on small system sizes. It is found that, unlike the previous studied 
case of equal inter and intra-chain coupling constants, the XY and Ising-like 
excitations of the quantum model decouple for large interchain coupling, giv- 
ing rise to pure Ising model critical behavior for the chirality order parameter 
and a superconductor-insulator transition in the universality class of the 2D 
classical XY model. 
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The one-dimensional frustrated quantum XY model (ID FQXY) has been introduced as a 
model for studying charging effects in a ladder of Josephson junctions in a magnetic field cor- 
responding to half a flux quantum per unit celjil. These charging effects arise from the small 
capacitance of the grains making up the ladder and leads to strong quantum fluctuations of 
the phase of the superconducting order parameter. As a result of the competition between 
the charging energy and the Josephson coupling between the superconducting grains, the 
one-dimensional array undergoes a superconductor to insulator transition at zero tempera- 
ture for decreasing capacitanceHi'S. The universality class of this transition is currently a 
problem of great interest specially in relation to experiments on two-dimensional supercon- 
ducting fllms and Josephson-j unction arraysi'ii'i0. For a chain of Josephson junctions this 
critical behavior has been identifled with that of the well known classical two-dimensional 
XY model!. However, while for the case of two coupled chains forming a ladder, in the 
absence of a magnetic fleld, this critical behavior can be shown to remain unchanged, in the 
presence of a magnetic fleld the behavior is more complicatedBI^. In particular, at half flux 
quantum per plaquette, corresponding to the ID FQXY mo del0 we study here, the critical 
beha vior is expected to be described by the 2D classical XY-Ising mo The existence 

of both XY and Ising-like excitations in this case result from the frustration induced by the 
magnetic fleld and are associated with the continuous U{1) symmetry of the phases of the 
superconducting order parameter and the plaquette chiralities which measures the direction 
of circulating currents in the Josephson-j unction ladder. 

The ID FQXY model is deflned by the Hamiltonian@ 



and consists of a one-dimensional chain of frustrated plaquettes as indicated in Fig. 1. 
The flrst term in Eq. (1) describes quantum fluctuations induced by the charging energy 
Ec = 4e^/C of a non-neutral superconducting grain located at site r, where e is the elec- 
tronic charge and C is the effective capacitance of the grain. The second term is the usual 
Josephson-j unction coupling between nearest-neighbor grains. 6^ represents the phase of the 
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superconducting order parameter and the couphngs Err' satisfy the Villain's 'odd rule' in 
which the number of negative bonds in an elementary cell is oddEl. This rule is a direct 
consequence of the constraint that, for the half flux case, the line integral of the vector po- 
tential due to the applied magnetic field should be equal to vr in units of the flux quantumi. 
In the classical limit {E^ = 0), the ground state of Eq. (1) has a discrete Z2 symmetry 
associated with an antiferromagnetic pattern of plaquette chiralities Xp = measuring 
the two opposite directions of the super-current circulating in each plaquette. In a previous 
work0, a particular case of the ID FQXY model, i.e., Ey = E^, has been studied in some 
detail and it has been found that in fact its critical behavior is consistent with the results for 
the 2D classical XY-Ising mo del§. From the critical exponents associated to the chirality 
order parameter the critical behavior has been identified as the one along the line of single 
transitions where both phase coherence and chiral order are lost simultaneously. However, 
the XY-Ising model has in addition to this transition line, two other branches correspond- 
ing to separate XY and Ising critical behavior which join the line of single transition at a 
bifurcation point located at some place in the phase diagram. The ID FQXY model studied 
previously corresponds to a particular path through this phase diagram; the one located in 
the region of single transitions. 

In this work we consider a generalized version of the ID FQXY model in which the 
inter (E^) and intra-chain (Ey) couplings constants have different strengths. In terms of a 
fluctuation-effective action which can obtained from an imaginary-time path integral rep- 
resentation of Eq. (1), the ratio between the couplings constants E^/Ey can be used to 
tune the system through the bifurcation point in the XY-Ising mo dei. In particular, for 
Ere » Ey the ID FQXY model is expected to have two separate transitions. In this work 
we find that in fact for E^/Ey ~ 2 the single transition found previously does decouple 
into two separated transitions. Using a Monte Carlo transfer matrix technique^ applied to 
the path-integral representation of the model we study the critical behavior of the chirality 
order parameter at a particular value of this ratio, E^/Ey = 3. We find, from a finite-size 
scaling analysis of extensive calculations on small system sizes, that the critical exponents 



are consistent with pure Ising model critical behavior as expected from the results for the 
XY-Ising model. Thus the superconductor to insulator transition in the related Josephson 
junction ladder is in the universality class of the 2D classical XY model. 

To study the critical behavior of the ID FQXY model, we find it convenient to use 
an imaginary-time path-integral formulation of the mo del0. In this formulation, the one- 
dimensional quantum problem maps into a 2D classical statistical mechanics problem where 
the ground state energy of the quantum model of finite size L corresponds to the reduced free 
energy per unit length of the classical model defined on an infinite strip of width L along the 
imaginary time direction, where the time axis r is discretized in slices Ar. After scaling the 
time slices appropriately in order to get a space-time isotropic model, the resulting classical 
partition function is given by Z = tre^^ where the reduced classical Hamiltonian is defined 
as 

H = —a ^ [ COS^Ot-j — OtJ+i) + COs(6't- j — ^^r+lj) 
r,j 

- COs{(f)r,j - (pT,j+l) + COs{(j)rJ - 0r+lj) 

+ t C0S(^.,, - <Pr,j) ] (2) 

In the above equation, 6 and denote the phases on the left and right columns in Fig. 1, 
and a = (Ey/EcY^'^ plays the role of an inverse temperature in the 2D classical model. 

One can now carry out a detailed study of the scaling behavior of the energy gap for kink 
excitations (chiral domain walls) of the ID FQXY model by noting that this corresponds 
to the interface free energy of an infinite strip in the model of Eq. (2). For large a (small 
charging energy Ec), there is a gap for creation of kinks in the antiferromagnetic pattern of 
Xp and the ground state has long-range chiral order. At some critical value of a, chiral order 
is destroyed by kink excitations, with an energy gap vanishing as |a — ad'^, which defines the 
correlation length exponent u. Right at this critical point, the correlation function decays as 
a power law < XpXp' >= \p ~ p'\~^ with a critical exponent rj. However to proceed further, 
one has to be able to calculate the free energy per unit length f{a) of the Hamiltonian on 
the infinite strip, which is usually obtained from the largest eigenvalue Ao of the transfer 
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matrix between different time slices as / = — In Aq. Here, the major difficulty in performing 
this type of calculation comes from the continuous degrees of freedom of the model which 
prevents an exact diagonalization of the transfer matrix. This can be overcome by using a 
Monte Carlo transfer-matrix method^ which has been shown to lead to accurate estimates 
of the largest eigenvalue even for this type of problems. Here we just summarize the main 
steps. The method is a stochastic implementation of the well-known power method to obtain 
the dominant eigenvalue of a matrix. First helical boundary conditions are implemented, 
in order to get a sparse matrix. Then, a sequency of random walkers Ri, 1 < i < r, 
representing the configurations of a column with L spins in the infinite strip is introduced 
with corresponding weights Wi. The number of walkers r is maintained within a few percent 
of a target value by adjusting the weights properly. A matrix multiplication can be 
regarded as a transition process with a probability density defined from the elements of 
the transfer matrix. In the procedure, a Monte Carlo (MC) step consists of a complete 
sweep over all random walkers and after a large number of MC steps an estimate of the 
largest eigenvalue can be obtained from the ratio between the total weights J2i of two 
successive MC steps. The implementation and some of the difficulties of the method are 
similar to the case of the two-dimensional frustrated classical XY mo dellli and the reader 
should refer to that work for further details. For the calculations discussed in this work, 
typically = 20000 random walkers and 80000 MC steps were used which correspond to 
1.5 X 10® attempts per (6', 0) pair. 

The interfacial energy for domain walls in the model of Eq. (2) can be obtained from the 
differences between the free energies for the infinite strip with and without a wall. However, 
because of the antiferromagnetic pattern of the chiralities Xp = only strips with an 
odd number of sites L will have a domain wall. Since one is required to obtain the free 
energy differences at the same value of L, we need to resort to an interpolation scheme for 
successive odd or even L to determine the interfacial free energy. Results for the interfacial 
free energy, defined as AF{a, L) = L'^Af{a, L), near the transition point etc, for 6 < L < 14, 
are indicated in Fig. 2 for a particular value of the ratio E^/Ey = 3. As in the previous 
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wor to obtain the critical exponents and critical temperature we now employ the finite-size 
scaling 

AF{a,L) = AiL'/'Sa) (3) 

where A is a scaling function and 6 = a — ac- In a linear ap proximation for the argument 
of A, we have 

AF{a,L) = a + bL^/''6a (4) 

which can be used to determine the critical coupling etc and the exponent u independently. 
The change from an increasing trend with L to a decreasing trend provides and estimate of 
ac, which from Fig. 2 gives ac = 1.16(2). Once the critical coupling is known, the correlation 
function exponent rj can obtained from the universal amplitude a in Eq. (4) through a result 
from conformal invariancelll, a = irrj, from which we estimate t] = 0.27(3). To estimate the 
correlation length exponent u we first obtain the derivative 5* = OAF /da near Oc, then it 
can easily be seeing that a log-log plot of S vs L gives an estimate oi 1/v without requiring 
a precise determination of a^- Of course, this is only valid for the linear approximation of 
Eq. (4). Assuming the data is in fact in this regime we obtain the result for S as indicated 
in Fig. 3 and get the estimate v = 1.05(6). The results for the critical exponents u and 
rj are good agreement with pure 2D Ising values u = 1 and 1] = 0.25 indicating pure Ising 
behavior. Moreover, this implies from the relation between the ID FQXY model and the 
2D classical XY-Ising|lll that the XY and Ising-like excitations have decoupled in this region. 

To show that in fact one has two decoupled and at the time separated transitions we 
now consider the results for the helicity modulus which measures the response of the system 
to an imposed phase twist. In the incoherent phase this quantity should vanishes while it 
should be finite in the coherent phase. The helicity modulus is related to the free-energy 
differences AF between strips with and without and additional phase mismatch of tt along 
the strip and is given by 7 = 2AF/7r^ for large system sizes. If the model is decoupled 
then the coherent to incoherent (or superconductor-insulator) transition should be in the 
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universality class of the 2D XY model, where one knows that the transition is associated 
with a universal jump of 27r in the helicity modulusEl. Finite-size effects smooth out this 
behavior as indicated in Fig 4 but the critical coupling can be estimated as the value of a 
at which AF = vr. This criteria leads to the estimate ac = 1.29 which is to be compared 
with the critical coupling for the destruction of chiral order found above, ttc = 1.16. This 
clearly indicates the transitions are well separated and thus one expects they are decoupled. 
In contrast, for the case of equal inter and intra-chain coupling constants studied previously 
these estimates were found to be in fact consistent with each otheiE. We have also performed 
less detailed calculations at other values of the ratio E^/Ey from which we can estimate that 
the Ising and XY transition merge into a single transition roughly at E^/Ey ~ 2. Since, the 
superconductor to insulator transition is to be identified with the loss of phase coherence! 
we reach the interesting result that in the ID FQXY, or alternatively, a Josephson-junction 
ladder, the universality class of the superconductor-insulator transition depends on the ratio 
between inter and intra-chain couplings. 

In conclusion, we have studied a generalized version of the one- dimensional frustrated 
quantum XY model which consisted in allowing for different strengths for the inter and 
intra-chain couplings constants. The model can be physically realized as a one-dimensional 
array of Josephson junctions in the form of a ladder and in the presence of an external 
magnetic field corresponding to a half fiux quantum per plaquette. It is found that, unlike 
the previous studied case of equal inter and intra-chain couplings, the XY and Ising-like 
excitations decouple giving rise to pure Ising behavior for chirality order parameter and a 
superconductor-insulator transition in the universality class of the XY model. Since these 
arrays can currently be fabricated in any desired geometry and with well-controlled param- 
eters it is hoped that these results will serve to motivate experiments in these systems. 
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FIGURES 

FIG. 1. Schematic representation of the one-dimensional frustrated quantum XY model with 

inter (i?x) and intra-chain (i-Ey) coupling constants. The antifcrromagnetic ordering of chiralities 
Xp = ±1 is also indicated. 

FIG. 2. Finite-size scaling of the interfacial free energy AF(q!, L) = L^A/(a, L) for kink (chiral) 
excitations. 

FIG. 3. S = dAF{a,L)/da evaluated near the critical coupling ac- The slope of the straight 
line gives an estimate oi l/u. 

FIG. 4. Behavior of the interfacial free energy AF = L'^Af for a system of size L = 12 
resulting from an imposed phase twist of tt. Vertical arrows indicate the locations of the Ising 
and XY transitions and the horizontal arrow the value AF = tt from where the XY transition is 
located. The Ising transition is located from the finite-size scaling of the chiral order parameter 
(Fig. 2) as discussed in the text. 
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